Controlling Anderson localization in disordered photonic crystal waveguides 
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We prove Anderson localization in a disordered photonic crystal waveguide by measuring the 
ensemble-averaged localization length, ^, which is controlled by the dispersion of the photonic crystal 
waveguide. In such structures, (, shows a 10-fold variation between the fast- and the slow-light regime 
and, in the latter case, it becomes shorter than the sample length thus giving rise to strongly confined 
modes. The dispersive behavior of ^ demonstrates the close relation between Anderson localization 
and the photon density of states in disordered photonic crystals, which opens a promising route to 
controlling and exploiting Anderson localization for efficient light confinement. 

PACS numbers: (42.25.Dd, 42.25.Fx, 46.65.+g, 42.70.Qs) 
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Quantum optics and quantum information technolo- 
gies require enhancement of light-matter interaction by, 
for example, confining light in a highly engineered 
nanocavity Quite remarkably, an alternative route 
towards light confinement exploits multiple scattering of 
light in disordered photonic structures, as originally pro- 
posed for electrons by P. W. Anderson [2]. The mech- 
anism responsible for Anderson localization is wave in- 
terference and hence occurs not only for electrons but 
also for, e.g., microwaves [1], acoustic waves [1], and even 
Bose-Einstein condensated matter waves @ thus illumi- 
nating the multidisciplinarity of the research field. In the 
case of light, indications of three-dimensional (3D) An- 
derson localization have been observed in random dielec- 
tric materials like powders composed of particles with 
casual shapes and sizes [6]. In these systems no control 
can be exerted over the frequency or spatial extent of the 
localized modes. 

A promising proposal on how to control multiple scat- 
tering is to introduce disorder in photonic crystals Q: 
the interference of multiply scattered light is expected to 
form Anderson-localized modes that appear randomly in 
space but at frequencies in or near the photonic crystal 
band gap. The characteristic length scale of Anderson 
localization is the localization length ^, which is the ex- 
ponential decay length of the confined modes after av- 
eraging over many realizations of disorder. In one- and 
two-dimensional (ID, 2D) systems, localization occurs for 
any degree of disorder when the sample length exceeds 
^ [Sj, in which case the photonic conductor becomes an 
insulator. Contrary to non-dispersive systems, photonic 
crystals offer the possibility to alter the photon density of 
states (DOS) allowing to control macroscopic transport 
properties [9] or spontaneous emission of photons lOj . 
and in particular to achieve dispersive Anderson local- 
ization. 

In this Letter, we show how to control and accurately 
tune Anderson localization of light to a spectral range 
of particular interest using dispersion in photonic crys- 
tal. In particular, we use a ID disordered photonic crys- 



tal waveguide (PCW) to prove the close relation between 
Anderson localization of light and the DOS. Indeed, pho- 
tonic crystals owe their success to the fact that their DOS 
can be accurately tailored a priori and we show how this 
property can be exploited to tailor Anderson localiza- 
tion. By explicitly measuring the dispersive localization 
length, we unambiguously demonstrate that the strongly 
confined modes appearing in the slow-light regime of 
PCWs ll[ are due to ID Anderson localization, which 
has been questioned recently In addition, we inter- 
pret our experimental data on the wavelength-dependent 
localization length with a model for the DOS of a PCW 
thereby explicitly linking the localization length and the 
DOS. 

Our samples consist of a membrane with a high refrac- 
tive index material (GaAs, n = 3.54) in which light is 
confined by total internal refiection. An ordered lattice 
of holes is etched in the structure forming a 2D photonic 
band gap that suppresses in-plane propagation of light. 
A waveguide is engineered in the structure by leaving out 
a row of holes. Light propagation in an ideal PCW is de- 
scribed by Bloch modes with a dispersion relation a;(fc), 
where uj is the wave frequency and k is the wave vec- 
tor. PCWs can generate slow light, i.e., light with a very 
low group velocity, Wg = duj/dk = cjn^, where c is the 
speed of light in vacuum. The slow-down factor rig, also 
referred to as the group index, is directly proportional to 
the DOS per unit length of the propagating mode in the 
PCW: DOS = (l/7r)afc/aa; = njii^c). Imperfections in 
a PCW lead to multiple scattering of light and the dis- 
persion relation breaks down [13] thus inducing localized 
modes in the slow-light regime where the sample length 
exceeds ^. 

To investigate the impact of disorder on the light prop- 
agation in PCW we randomly vary the hole positions 
in the three rows above and below the waveguide us- 
ing a Gaussian random number generator function (Box- 
Muller). The degree of disorder in each sample, (5, is 
characterized by the standard deviation of the hole po- 
sition with respect to the lattice constant varying from 
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FIG. 1: (Color online) (a) and (b) Scanning electron micro- 
scope images of photonic crystal waveguides (top- view) with 
engineered disorder on the positions of the holes of (5 = 0% 
and S = 12% (standard deviation relative to the lattice con- 
stant), (c) Setup used to measure the localization length. 
The tip of a single mode tapered fiber (F) is placed on the 
photonic crystal waveguide. TE-polarized light scattered from 
the waveguide is collected with an optical microscope objec- 
tive (M), a spatial filter (SF), and a polarizer (P) and sent to 
a spectrometer (S). (d) Measurement of strongly fluctuating 
light intensity in the slow-light regime of a PCW with S = 0% 
versus the distance from the light source. The inset of the fig- 
ure shows a spectrum taken at L = 150 \im from the fiber 
tip. 



5 ^ 0% to 5 ^ 12% (Fig. ma), (b)). The samples con- 
sist of a triangular lattice of holes with a lattice con- 
stant a = 240 nm, a filling fraction / = 0.330 ± 0.006, 
membrane height h — 160 ± 5 nm, and a total length 
Lq = 1 mm. 

To characterize Anderson localization in PCWs we use 
the optical setup illustrated in Fig. [IJc). A continuous 
wave Ti:sapphire laser tuneable within A — 700 — 1000 nm 
is coupled into a single mode tapered fiber with a tip 
diameter comparable to the waveguide width. The fiber 
evanescent mode couples to the waveguide mode by plac- 
ing the fiber tip in close proximity of the PCW. We mea- 
sure light scattered out-of-plane from the PCW with a 
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FIG. 2: (Color online) (a) Exponential fit of the ensemble- 
averaged light transmission in PCWs with S = 0% as a func- 
tion of the distance from the source for two wavelengths: 
Ai = 890 nm (gray points) and A2 = 916 nm (black points). 
The spectral resolution is 2 nm and the slope equals the in- 
verse of the localization length, (b) The calculated density 
of states (DOS) and group index (rig) of an ideal structure 
without disorder. The two specific wavelengths Ai (gray line) 
and A2 (black line) are indicated that lie in the fast- and 
slow-light regimes, respectively. The calculated spectral posi- 
tion of the cutoff of the ideal waveguide mode gives the best 
fit to our experimental data (explained in the text), (c) Plot 
of the reduced chi-square, x^, for the exponential fit to the 
ensemble-averaged light transmission in the fast- (gray trian- 
gles) and slow-light regime (black squares) versus number of 
realizations of disorder. The dashed lines are guides to the 
eyes. 



high amplification microscope objective (NA — 0.8) as a 
function of the wavelength and the distance L from the 
light source, i.e., the fiber tip. The measurement starts at 
L = 150 |J.m from the fiber tip to avoid any spurious effect 
due to the evanescent field or light that is not coupled to 
the waveguide mode. Fig. [Ijd) shows a measurement of 
scattered light intensity versus L in the slow-light regime 
(at A = 916 nm where Ug > 15) for a single realization 
of disorder. In this case, S = 0%, the PCW is only af- 
fected by intrinsic unavoidable disorder introduced in the 
fabrication process. The inset of the figure shows a spec- 
trum of the scattered light intensity at L = 150 |^m. The 
strong fiuctuation in the light intensity is a signature of 
ID Anderson-localized modes. The modes appear to be 
spatially and spectrally separated, which constitutes a 
criterion for Anderson localization [l3 |. 

In the Anderson localization regime the fluctuat- 
ing light transmission, T{L), decays exponentially after 
ensemble-averaging, giving ln(T) — — TE*]. We have 
used this property to extract the localization length and 
confirm the Anderson localization criterion of ^ < L. To 
do this, we perform ensemble averaging by measuring the 
light leakage versus L and wavelength for 8 different re- 
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FIG. 3: (Color online) Plot of the localization length, ^, of 
the PCW with S — 0% as a function of the wavelength and the 
group index of the ideal structure ng. A strongly dispersive 
localization length is found that reaches a minimum in the 
slow-light regime (rig > 15). The scaling ^ oc n"^ (red curve) 
is obtained from a model of the light scattering cross section 
in the PCW. 



alizations of disorder. In particular, we probe different 
spatial realizations of disorder by moving the fiber to- 
gether with the microscope objective along the waveg- 
uide and repeating the measurement. Fig. [2ja) shows 
the linear fit of ln(T) in the sample with 6 = 0% in 
the fast-light regime ?^g(Ai) < 15 and slow-light regime 
ng(A2) > 15 (corresponding to low and high DOS, re- 
spectively, see Fig. [5Jb)). For these spectral positions we 
extract the localization lengths C(Ai) = (263 ± 28) |^m 
and ^(A2) = (27 ± 2) p.m. We obtain rig from numeri- 
cal simulations of the disorder-free ideal structure, which 
quantifies the spectral position with respect to the cut- 
off of the ideal waveguide mode. This spectral position is 
obtained from the best fit to our experimental data with 
a ID model (explained later on). It is within the uncer- 
tainty in the waveguide mode band edge position, which 
is caused by uncertainties in the fabrication parameters 
given previously. 

Fig.lHc) shows the goodness of the fits, x^, performed 
in Fig. [2Ja) to test the degree of convergence of the en- 
semble average to a single-exponential decay. The fluctu- 
ations in the data are due to speckles not fully ensemble- 
averaged. A disorder realization in our experiment con- 
sists of a scan of T{L) varying L over 130 |J.m. Thus, the 
number of different disorder realizations that we can per- 
form without adding repeated statistics to the ensemble 
average is limited by the sample length, Lq = I mm. Fur- 
thermore, the goodness of the fit does not reach the opti- 
mum value = 1 corresponding to a single-exponential 



decay, which may be due to different reasons. Deviations 
from perfect verticality of the air holes can break the 
symmetry in the out-of-plane direction and couple TE- to 
TM-polarized waveguide modes giving rise to strongly lo- 
calized modes [l^ . This polarization mixing mechanism 
may lead to multi-exponential decay of the ensemble- 
averaged transmission that is not resolved in the present 
experiment. 

Fig. [3] contains the main contribution of this Letter: it 
shows the dispersive behavior of the measured localiza- 
tion length. We observe a 10-fold variation in the localiza- 
tion length with wavelength. In the slow-light regime, ^ 
reaches its minimum value of ~ (27 ± 2) |j.m and becomes 
much smaller than Lq, giving rise to strongly Anderson- 
localized modes. The dispersive behavior of the localiza- 
tion length gives directly control over the extension of 
the modes, which can be precisely tuned by varying the 
wavelength. We can also tune the waveguide mode by 
varying the fabrication parameters (typically a and /) 
thus, controlling the frequency of localized modes. 

Any loss mechanisms of the light trapped in the PCW 
will in general influence the measured localization length 
depending on the spatial length scale of the loss process. 
In presence of light losses, we have ln(T) = —L/2£ where 
the decay length i is defined as = + £i^^ + £e^^, 
where £i is the material inelastic absorption length and 
£e is the extinction length associated with out-of-plane 
losses. GaAs has a very low optical absorption coefficient 
(< 100 cm~^ at a wavelength of A = 915 nm) correspond- 
ing to £i > 1 m. This value might be reduced by surface 
effects at the holes of the photonic crystal but is still 
expected to be much larger than Lq. Furthermore, we 
quantify £e applying our fabrication parameters to re- 
cent 3D numerical simulations of Bloch-mode scattering 
in PCWs [13] ■ In particular, for 6 = 0%, for which the 
standard deviation of the hole positions is smaller than 
2 nm, we obtain ig 400 mm. Both £i and £e are much 
larger than Lq, thus not affecting the localization length 
we extract from our data, i.e., we can approximate £ ^ ^ 
at least for weak disorder. 

In the following, we model the function £(ng). The re- 
lation N ■ £s holds for ID systems la, [l8|, where N 
is the number of electromagnetic modes that the system 
can sustain and is the scattering mean free path. Our 
samples are single-mode PCWs and hence the localiza- 
tion length is ^ ~ 4- The dispersive behavior of 4 in 
3D photonic crystals has been measured and explained 
in terms of the DOS very recently [l^. Two separate 
mechanisms determine the excitation of the scatterer 
and the subsequent radiation from the scatterer [2^ . The 
couphng to the scatterer is described by the DOS of the 
mode of the excitation beam [2l|| , which in our system is 
proportional to the group index of the waveguide mode. 
The second process is described by the local density of 
states (LDOS). Ignoring the minor contributions of cou- 
pling to unconfined radiation modes, the LDOS is also 
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determined by the group index of the waveguide 2^, l23| . 
This apphes to every scattering event giving rise to a 
modified scattering cross section a in PCWs scahng as 
a cx ng^(a;). The scattering mean free path in a random 
medium can be expressed as 4 = l/ftc [2^], where ps 
is the density of scatterers. For ID single-mode PCWs 
we therefore predict ^ oc ng~^(a;), which is in very good 
agreement with our experimental data in the fast-light 
regime. The red curve in Fig. [3] represents the best fit 
to our experimental results. It gives the spectral position 
of the waveguide mode cutoff used in Fig. ^h). The fit- 
ting parameters are the number of data fitted and the 
spectral position of the cutoff of the waveguide mode. 
The same scaling of ^ oc ng~^ can also be recovered from 
ID random matrix theory [2^ and it is confirmed by 3D 
numerical simulations in PCWs . Our model explains 
not only the dispersion in the localization length but also 
the scaling of light losses with group index in PCWs as 
Ug~^ (ng^). Such a scaling has been systematically ob- 
served in the literature [26| and our model provides a 
physical explanation. From the data in Fig. [3] we observe 
that the model breaks down deep in the slow-light regime 
where strongly localized Anderson modes appear. This is 
expected since the model is based on the calculated DOS 
of the ideal structure without disorder, which is modified 
in the regime of strong multiple scattering and Anderson 
localization. 

Finally, Fig. 0] plots the measurement of the localiza- 
tion length in samples with increasing amount of disor- 
der and the mean value of the localization length, ^. The 
latter decreases with disorder, reaches a minimum for 
S = 6%, and increases for 6 > 6%. This behavior is a clear 
proof that losses are not dominant in our experiment even 
for strong disorder where £e is predicted to shorten 17j . 
The increase of ^ for (5 > 6% cannot be explained by 
an increase of light losses (decrease of ie)- That would 
lead to a decrease of the measured decay length as op- 
posed to our observations. The behavior of ^ for weak 
disorder is predicted in Ref. 
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17| and we propose here a 



possible explanation for the increase observed for large 
amounts of disorder. Thus, the local disorder introduced 
in our samples only in the three rows above and below 
the waveguide could imply that several spatial modes, 
iV, are effectively introduced when increasing disorder. 
This would increase the localization length according to 
^ oc N -is effectively decoupling ^ and £s and therefore 
enable a regime of light diffusion where £s < L < S^. 
This could open a new route to investigate the crossover 
between diffusion and localization regimes in a disper- 
sive quasi-lD system, which is not possible in standard 
quasi- ID disordered systems isl 27|. 

In conclusion, we have demonstrated the close rela- 
tion between Anderson localization of light and the elec- 
tromagnetic DOS in PCWs. We explain multiple light 
scattering in PCWs with a DOS dependent scattering 
cross section explaining the strong dispersion of the lo- 
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FIG. 4: (Color online) Localization length, ^, as a function of 
the group index, ng, measured in PCWs with various amounts 
of disorder 5 = 0% (o), 3% (□), 6% (A), 9% (*), and 12% 
(o). The inset shows the mean value of ^ as a function of 5. 



calization length and the appearance of strongly con- 
fined modes in the slow-light regime. These results are 
of fundamental importance since they impose limitations 
for slow-light devices based on PCWs such as single- 
photon sources At the contrary the strongly con- 
fined Anderson-localized cavities with tailored properties 
appear very appealing candidates for experiments on cav- 
ity quantum electrodynamics j28] or random lasing. 
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